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Abstract 

Monte Carlo simulation is one of the main applications involving the use of random 
number generators. It is also one of the best methods of testing the randomness prop- 
erties of such generators, by comparing results of simulations using different generators 
with each other, or with analytic results. Here we compare the performance of some 
popular random number generators by high precision Monte Carlo simulation of the 2-d 
Ising model, for which exact results are known, using the Metropolis, Swendsen-Wang, 
and Wolff Monte Carlo algorithms. Many widely used generators that perform well in 
standard statistical tests are shown to fail these Monte Carlo tests. 



1 Introduction 



Monte Carlo simulation is an important numerical technique for studying a wide range 
of problems in the physical sciences.^ Being a probabilistic technique, it relies heavily on 
the use of pseudo-random number generators-^' S @ The generation of random numbers on 
a computer is a notoriously difficult problem. An ideal random number generator would 
provide numbers that are uniformly distributed, uncorrelated, satisfy any statistical test 
of randomness, have a large period of repetition, can be changed by adjusting an initial 
"seed" value, are repeatable, portable, and can be generated rapidly using minimal computer 
memory. 

Many statistical tests have been developed to check for randomness^ ! and in most 
cases the period of the generator can be calculated (at least approximately). As noted in 
a number of recent review ar tides,!' & i random number generators provided by computer 
vendors or recommended in computer science texts often have been (and unfortunately 
continue to be) of poor quality. Even generators that perform well in standard statistical 
tests for randomness may be unreliable in certain applications, as has been found in some 
Monte Carlo simulations B I !' 0' @ 

There are two main types of random number generators for producing sequences of 
pseudo-random integers Xf. 

1. Linear congruential generators (LCGsjl' i 

Xi = A* Xj_i + B mod M 

which we will denote by h(A,B,M). The period is M for suitably chosen A and B 
(M- 1 if B = 0). 

2. Lagged Fibonacci generators (LFGs)& ^ 

Xi = Xi_p Xi-Q 

which we will denote by F(P, Q, 0), P > Q, where is any binary arithmetic oper- 
ation, such as +, — , * or © (the bitwise exclusive OR function XOR). The arithmetic 
operations are done modulo any large integer value, or modulo 1 if the X's are repre- 
sented as floating point numbers in the interval [0,1), as can be done if the operation 
is + or — . Multiplication is done on the set of odd integers. For 6-bit precision X's, 
the period is (2 P — l)2 b ^ 1 , or (2 P — 1)2 & ~ 3 for multiplication, for suitably chosen lags.H 
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It is possible to find sets of parameters (A,B,M) or (P,Q,Q) for which these two 
types of generators work well for most practical purposes, and it is possible to improve the 
performance of these generators by increasing M or P.i There are practical limits on these 
two parameters: M should not be very much greater than machine precision to avoid using 
slow multi-precision arithmetic, and a large lag P means storing a large array of previous 
numbers in the sequence (the "lag table") which may be subject to memory constraints. 
However on most modern computers adequate values of M and P can be found which are 
well within these limits. 

Linear congruential generators have two major defects. The first is that the least sig- 
nificant bits of the numbers produced are highly correlated, and a resultant "scatter-plot" 
of ordered pairs of random floating point numbers in the interval [0,1) shows regular lattice 
structure.^ 00' lil They are also known to have long-range correlations, especially for 
intervals which are a power 

f 2 JB 1 0. § Another problem is that for 32-bit integers the 
period of these generators is at most 2 32 , or of order 10 9 . On a modern RISC workstation 
capable of around 10 8 floating point operations per second, this period can be exhausted 
in a matter of minutes. This can be alleviated by the use of 64-bit precision, however the 
correlation problems still remain (although to a lesser degree). In spite of these problems, 
LCGs with well-chosen parameters perform well in most standard statistical tests, and an 
LCG (unfortunately not always with well-chosen parameters!) is provided as the default 
generator on many computer systems. 

Lagged Fibonacci generators using arithmetic operations (+,—,*) give good results 
in standard statistical tests with very modest lags on the order of tens.i When the binary 
operation used is XOR, these generators are referred to as generalized feedback shift register 
generators .0' Marsaglia has shown that XOR is one of the worst operations one can 
use in a generator of this type, and strongly recommends the use of standard arithmetic 
operations that have much longer periods and perform much better in statistical tests.^ 
Although shift register generators pass statistical tests when the lag is large enough (of 
order hundreds), & ^ 1^1 very little (apart from the period) is known theoretically about 
these generators, and they have produced biased results in Monte Carlo studies of the Ising 
model in two0 and three§ dimensions, and of self-avoiding random walks.0' ^ 

Mixing two different generators is believed to improve performance in some cases,^' ^ 
and many generators that perform well in statistical tests are of this kind. Marsaglia has 
suggested a fast, simple Weyl (or arithmetic sequence) generator! — H — I 

Xi = - K mod M, 



2 



with K a constant relatively prime to M, that can be effectively combined with a lagged 
Fibonacci generator. Adding a Weyl generator also increases the period of the combined 
generator by a factor of M (the period of the Weyl generator). L'Ecuyer^l has shown how 
to combine two different 32-bit LCGs to produce a mixed generator that passes the scatter- 
plot test and has a long period of around 10 18 , thus overcoming some of the drawbacks of 
standard LCGs. Although these mixed generators perform well in empirical tests, there is 
little theoretical understanding of their behavior, and it is quite possible that mixing two 
generators may introduce new defects of which we are unaware. A good single generator 
may therefore be preferable to a mixed generator. 

LCGs have the advantage that we have a relatively good (although still limited) theoret- 
ically understanding of their randomness properties. They are known to be defective, but 
their defects are fairly well understood (for example, the lattice structure of an LCG can 
be determined analytically using the spectral testP), and in practice they work quite well. 
There is clearly a need for better random number generators, and LFGs and mixed genera- 
tors are prime candidates. However currently there is little or no theoretical understanding 
of these and other generators, and they are used mainly on the basis of their performance 
in statistical tests. They are believed to overcome some of the flaws of LCGs, although 
this has not been proven and they may possess other flaws of which we are unaware. It is 
therefore extremely important to subject random number generators to a wide variety of 
precise statistical tests. 

2 Monte Carlo Tests 

One practical way to test a random number generator is to use it for Monte Carlo simulation 
of the two dimensional Ising model.0 This simple model has been solved exactly for a finite 
latticej^l so that values of the energy and the specific heat (the variance of the energy) of 
the system calculated from the Monte Carlo simulation can be compared with the known 
exact values. 

A number of different Monte Carlo algorithms can be used to simulate the Ising model. 
Here we will concentrate on the three most widely used methods: the Metropolis algorithm,0 lil 
which updates a single site of the lattice at a time; the Swendsen-Wang algorithm,^ which 
forms clusters of sites to be updated collectively; and the Wolff algorithm^ which updates 
a single cluster of sites. Each of these algorithms uses random numbers in a very different 
way. The Swendsen-Wang and Wolff cluster update algorithms are extremely efficient and 
allow very precise Monte Carlo simulations of the Ising model, easily reducing statistical 
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errors in the energy to better than one part in 10 5 . This precision provides us with a very 
effective practical test of the randomness of a pseudo-random number generator, and in 
particular its suitability for Monte Carlo simulation. 

Ferrenberg et recently showed that some "good" random number generators, which 
perform well in standard statistical tests, fail the "Monte Carlo test" ; that is, they produce 
incorrect results when used in Monte Carlo simulations of the Ising model, especially using 
the Wolff algorithm. The generators studied by Ferrenberg et al. were: 

i. CONG, the linear congruential generator L(16807, 0, 2 31 - 1).§ I 

ii. Two shift register generators, F(250, 103, ©) and F(1279, 1063, ©)S 

iii. SWC, a subtract-with-carry generator based on F(43, 22, — ).^ 

iv. SWCW, a combined subtract-with-carry and Weyl generator .0 

In spite of the premise of that paper, CONG and the shift register generators are in fact 
known to be not good random number generators. CONG has been recommended by a 
number of authord!^ as one of the best 32-bit linear congruential generators, however it 
still suffers the small period and correlated low order bits of these generators. Shift register 
generators have been criticized by Marsaglia, who showed that those with small lags (less 
than 100) performed poorly in statistical tests.il However similar tests of F(250, 103, ©) gave 
good results]^ 11 and Kirkpatrick and Stoll also obtained reasonable results with Monte 
Carlo testsS 

Subtract-with-carry generators are another variation of LFGs, where the standard op- 
eration of subtraction is replaced by subtraction with a carry bit C, as follows: 

Xi = Xi-p — Xi^Q — C, 

if Xi > 0, C = 0, 

if Xi < 0, Xi = Xi + M, C = 1. 

This greatly increases the period of the LFG, to M p — for suitably chosen P, Q and 
M,0 compared to approximately M2 P for a comparable LFG using subtraction. We have 
used M = 2 32 — 5, which gives very long periods for modest lags. Although advocated 
by Marsaglia,^ there were no known published results on statistical tests of the SWC or 
SWCW generators prior to the results of Ferrenberg et al. , so again there was little support 
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for their claim that these are "good" generators. Recently the shift-with-carry generators 
were in fact shown to perform poorly in standard statistical tests @ 

In this paper the work of Ferrenberg et al. has been extended by studies of both the 
"good" generators of that paper, and some "better" generators, which are listed below. In 
this work there are also more, and in some cases longer, independent runs for each generator, 
to obtain better error estimates and to better explore the effect of different initial seeds. 

In a recent review of random number generators,i James recommends 3 mixed genera- 
tors: 

1. RANECU, L'Ecuyer's mixed LCG combining L(40014,0, 2147483563) and 
L(40692,0, 2147483399)0 

2. RANMAR, Marsaglia's combined LFG F(97,33, -) and Weyl generator.il 

3. RCARRY, a subtract-with-carry generator^ based on F(24, 10, — ) (this is the same 
as SWC but with a smaller lag). 

We also tested the above generators, plus the following: 

4. RAND, the default 32-bit C and Unix generator L(1103515245, 12345, 2 31 - 1). 

5. DRAND48, another standard C and Unix generator with larger modulus and period, 
based on L(5DEECE66D W , B m ,2 i8 ). 

6. RANF, another 48-bit LCG, L(2875A2E7B175 W , 0, 2 48 ), which is the standard gen- 
erator used on CRAY and CDC CYBER machines.il 

7. RAN2, which is RANECU augmented by shuffling the order of the output values. ^1 

8. LFGs of different lags, using +, — , * and ©. 

9. LFGs using + and © with 4 "taps" S S ^ i.e. 

Xi = Xi_p Aj_Q Xi_R Aj_s, 

which we will denote by F(P, Q, R, S, 0). 

We followed Marsaglia and James by initializing each bit of the seed tables in the LFGs 
by using a combination LFG and LCG (see the routines RSTART in Ref. 28 and RMARIN 
in Ref. 3). We also tried using RAND to initialize every element of the seed tables, or every 
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bit of every element in the seed tables, which had little or no effect on the quality of the 
LFGs. 

For each random number generator, 25 independent simulation runs with different initial 
seeds were performed, on a network of IBM RS/6000, HP Apollo 9000, and DEC 5000 
workstations. Each simulation was between 10 6 and 5 x 10 7 sweeps of a 16 x 16 lattice at 
the critical point of the 2-d Ising modelS' The number of random numbers generated 
per sweep per site varies with the Monte Carlo algorithm used, with an average of 0.87 
for Metropolis, 0.93 for Wolff, and 1.85 for Swendsen-Wang. For the Metropolis algorithm 
we chose to visit the sites to be updated in order, rather than randomly, to provide a 
more effective way of probing any regularity or lattice structure in the sequence of random 
numbers, especially for the linear congruential style generators which are known to suffer 
from this problem! B 111 0> @- i 

Error estimates for each simulation were obtained by standard methods of binning the 
data, with a bin size much greater than the autocorrelation time.0 The error in the mean 
of the 25 combined results was also calculated, treating them as independent data sets. 
Two measures were used to compare the Monte Carlo results with the exact results: the 
deviation A between the mean of the combined results and the exact value as a multiple 
of the error in the mean a, and the chi squared per degree of freedom y 2 for the 25 data 
setsB The htst test check, tot an y h.as in the average over an rnns. the second checks 
for discrepancies in the statistical fluctuations expected between the individual runs. A 
generator is judged to have failed the Monte Carlo test if A > 3.3a, \ 2 > 2.0, or \ 2 < 0.34, 
all of which should occur with probability less than 0.001 for a truly random generator.il 

3 Results 

The results for A/a, the difference between the exact and simulated values of the energy 
and specific heat given as a multiple of the errors in the mean, are presented in Tables |l] and 
||, along with the values of x 2 - Failure of a test is indicated in bold type. The generators 
are grouped into 4 categories, determined by a different level of precision of the simulations. 
Table |Tj shows generators which we would classify as bad or very bad (at least for this type 
of Monte Carlo application). The very bad generators failed at least one of the tests with 
10 6 Monte Carlo sweeps per run, with the bad generators failing after 10 7 sweeps per run. 

Table |2| shows generators which we would classify as good or very good. The good 
generators failed one of the tests at a level of 5 x 10 7 sweeps for the Wolff and Metropolis 
algorithms, and 10 7 sweeps for the SW algorithm (which uses about twice as many random 
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numbers per sweep). The very good generators passed all the tests at this level, which 
involves generating on the order of 10 10 random numbers for each of the 25 independent 
simulations, or approximately 3 x 10 11 random numbers in total. In contrast, the errors 
caused by using very bad generators were generally apparent after using less than 10 9 
random numbers, in simulations which took only about an hour on a workstation. 

Fig. [l| shows the relative error in the specific heat for the Wolff algorithm versus the 
lag of the Fibonacci generator, for the binary operations addition, subtraction, and XOR. 
In all cases the XOR operation was about an order of magnitude worse than addition and 
subtraction. Since in Monte Carlo simulation an order of magnitude decrease in the error 
requires 100 times as many iterations, the difference between the quality of the LFG with 
different operations is substantial. Quite large lags of at least 1000 are required to reduce 
the error to less than 0.1%, however the percentage error for a given lag P goes roughly as 
e~ p ', so performance can be greatly improved with a moderate increase in the lag. For a 
lag of 4423 the generators gave correct results for all binary operations within the errors of 
the simulations. 

Table || compares the results for the Wolff algorithm for various generators based on 
F(43, 22, 0), where the binary operation is XOR, subtraction, subtract-with-carry, and mul- 
tiplication. The results of combining this lagged Fibonacci generator with a Weyl generator 
(as in SWCW or RANMAR) are also shown. We can see that the shift register generator 
using XOR performs very poorly, with errors of nearly 10% in the specific heat. Using 
subtraction performs an order of magnitude better, however adding a carry bit does not 
provide any extra improvement. Mixing in the Weyl generator reduces the errors by nearly 
another order of magnitude. Using multiplication instead of subtraction produces the most 
dramatic improvement, for little extra computational cost on modern RISC workstations. 
In Table |I| the standard 2-tap LFG is compared to a 4-tap version of the same lag, which 
gives substantially better results, as was seen by Ziff for self-avoiding random walks. 1 — 1 

The two 32-bit LCG generators both gave consistent results at the level of 10 6 sweeps, 
for which the number of random numbers required for each simulation is less than the 
period of these generators. Both failed the tests at the level of 10 7 sweeps, which requires 
producing about as many random numbers as the period. This suggests that the failure 
is due to the short period of these generators rather than the lack of randomness. This is 
supported in the case of RAND by the fact that some of the \ 2 values in Table |l| are smaller 
than expected, i.e. the deviations from the exact value of all the independent runs are too 
small. This is probably due to the fact that each run exhausts the period, so that different 
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runs are using similar sequences of random numbers and are therefore correlated to some 
extent. 

The mixed LCG generators RANECU and RAN2 were among the best generators, 
although they were also the slowest. This good performance was rather unfortunate in the 
case of the RAN2 generator, since the authors of Numerical Recipes have guaranteed RAN2 
to produce "perfect" random numbers, with perfect defined as "we will pay $1000 to anyone 
who convinces us otherwise (by finding a statistical test that RAN2 fails in a non-trivial 
way, excluding the ordinary limitations of a machine's floating point representation). "il 

The subtract-with-carry generators RCARRY and SWC were among the worst of the 
generators tested, which agrees with the results of Refs. 11 and 22. With the notable 
exception of the version using multiplication, the lagged Fibonacci generators performed 
very poorly for lags under 100 (under 1000 for the case of ©), and non-random effects were 
measurable even for lags of over 1000. In contrast, standard statistical tests by Marsaglia 
gave good results for LFGs using subtraction, even for lags less than 100 (except for the 
"birthday spacings" test). ^ ^ Marsaglia found that LFGs using multiplication performed 
very well in statistical tests even for small lags, and this is also true for the Monte Carlo 
tests, where multiplication gave by far the best performance for a given lag. Generators 
based on LFGs performed worst for the Wolff algorithm, with some small lag generators also 
failing the test with the Metropolis algorithm. LCGs performed worst on the Metropolis 
algorithm. 

GrassbergeiO tested F(250, 103, 0) using Monte Carlo simulations of random walks, 
and conjectured that this generator has large correlations over long times which should 
only be seen for Ising model simulations using lattices larger than 16 2 . We have also done 
simulations on a 128 2 lattice to compare the corresponding errors. The statistical error in 
the mean energy is 



where Tj„t is the integrated autocorrelation timep t-3 Ch is the specific heat, and V is the 
lattice volume. Since Tint f° r the Wolff algorithm has been measured to be 2.6 for V = 16 2 
and 8.2 for V = 128 2 ,S and the specific heat increases by 1.69 for the larger lattice,^ the 
statistical error in the mean energy will be approximately the same for 8.5 x 10 4 sweeps of 
the 128 2 lattice as for 10 6 sweeps of the 16 2 lattice, which was indeed found to be the case 
in our simulations. Table || shows that the discrepancy in the average energy caused by the 
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random number generator is actually much smaller for the larger lattice size. Since for the 
specific heat the statistical error increases even more rapidly with increasing lattice volume, 
smaller lattices seem to be more effective for testing some random number generators using 
Monte Carlo simulations of the Ising model. Of course the inverse result is also true - some 
random number generators will perform better in Monte Carlo simulations on large lattices. 

4 Conclusions 

Lagged Fibonacci generators using the operations of addition, subtraction or XOR (exclusive 
OR) can give poor performance, especially for the Wolff algorithm, unless the lag is very 
large. Using addition or subtraction gives substantially better performance than the shift 
register generators using XOR. Using multiplication gives extremely good performance even 
for small lags. Adding a carry bit to an LFG using subtraction (the subtract-with-carry 
generators) gives no improvement in the performance of these generators, however adding 
a simple Weyl generator greatly improves the quality of the LFG. 

The multiplicative lagged Fibonacci generator F(P, Q, *) was one of the best generators 
we tested. This generator showed excellent randomness properties even for very small 
lags, with only a slightly greater computational cost than an LCG, or (on modern RISC 
processors) an LFG using addition or subtraction. A multiplicative LFG can be given 
an arbitrarily large period by simply increasing the lag. A lag of only 43 gives a period 
of order 10 21 for 32-bit integer arithmetic, and extremely good randomness properties. 
The only drawback of these generators is the lack of a solid theoretical understanding of 
their properties. More theoretical studies and experimental tests should be done on these 
generators, since they appear to be very promising candidates for a good general purpose 
random number generator. 

The 32-bit linear congruential generators perform well up to the point where their pe- 
riod is exhausted, with RAND seemingly better than CONG. The 48-bit LCGs such as 
DRAND48 gave excellent results, and have a large enough period (of order 10 14 ) for most 
current applications. LCGs using even larger integers, such as L(13 13 , 0, 2 59 ), show very 
g ood perforce in standard sta.is.ica! resell and have even longer periods. These 
longer period LCGs usually require multi-precision arithmetic and are therefore relatively 
slow (e.g. DRAND48 is 6 times slower than RAND on a DECstation 5000), however they 
should become more popular in the near future, when 64-bit microprocessors become com- 
monplace. Apart from an increased period, large M LCGs also have better spectral (lattice) 
properties, however the correlations inherent in LCGs are still present. Combining a good 
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LCG with another generator, such as an LFG or another LCG (as with RANECU and 
RAN2), may further reduce (or even eliminate) these correlations, however it is possible 
that this may introduce other unknown defects. Again, we are hampered by the lack of 
a good theoretical understanding of these algorithms. In general it is probably advisable 
to stick with a good large M LCG, which should work perfectly well for most applica- 
tions. However it is known that these generators can perform poorly on vector and parallel 
computers, where the power-of-2 correlations can be accentuated. 

Note that by the year 2000 supercomputers will have Teraflop (10 12 floating point op- 
erations per second) performance, and a Teraflop-year of computation (3 x 10 19 flops) will 
become realizable for such problems as Monte Carlo simulation of lattice QCD and con- 
densed matter physics.^ It is therefore likely that large scale Monte Carlo simulations only 
ten years from now will exhaust the period (of roughly 10 18 ) of 64-bit LCGs or mixed 32- 
bit LCGs. However a 96-bit or 128-bit LCG, or a mixed generator made up of two 64-bit 
LCGs (similar to the RANECU generator studied here), should have both the randomness 
properties and the extremely large period necessary for any application in the forseeable 
future. These multi-precision arithmetic and mixed LCG algorithms are the slowest of the 
algorithms tested here, however it should be noted that the speed of a random number 
generator is often irrelevant, since in most applications the amount of time spent generat- 
ing the random numbers is insignificant compared to the rest of the calculation. In most 
applications the quality of the random numbers is far more important than the speed with 
which they are generated. 

Mixed lagged Fibonacci generators such as RANMAR have extremely long periods (10 43 
for RANMAR), however for high precision work the generator F(97, 33, — ) on which RAN- 
MAR is based should be replaced by a longer lag generator with better randomness prop- 
erties, such as F(250, 103,-), F(607, 273, -), or F(1279, 1033, -). The extra memory re- 
quirement is negligible for current workstations and high performance computers, except 
perhaps for fine grained massively parallel machines with limited memory per processor. 
Mixed generators offer a greatly increased period, and empirical tests indicate that they 
can have better randomness properties than the single generators on which they are based. 
The mixed generators were among the best tested here, however they are not as theoret- 
ically well understood as single generators, so it is possible that unexpected correlations 
may occur. They should therefore be used with caution. 

Our theoretical understanding of random number generators is quite limited, and no 
amount of statistical testing can ever determine the quality of a generator. It is there- 
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fore prudent in any stochastic simulation to use at least two very different generators (for 
example, a good large M LCG, a multiplicative LFG, or a good mixed generator such 
as RANMAR or RANECU) and compare the results obtained with each, in order to be 
confident that the random number generator is not introducing a bias in the results. 

Finally, we should note that it is unfortunate that most of the poorly performing gener- 
ators tested here are recommended in many texts and are available by default to the unwary 
user on many computer systems. @' i It should be no more acceptable for a computing envi- 
ronment to have a default random number generator that is known to be bad, than to have 
an incorrect implementation of a standard mathematical function. Since faster computers 
and better algorithms are improving the precision of Monte Carlo and other stochastic sim- 
ulations at a rapid pace, it is important to continue to search for better random number 
generators with very long periods, and to make more precise and varied tests of these gen- 
erators. This is particularly true for high performance computers with vector or parallel 
architectures, where methods for generating independent random numbers in parallel are 
required.^' ^ 
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Table 1: Results of Monte Carlo simulations of the 2-d Ising model using different random 
number generators. The first line for each generator shows the deviation of the Monte Carlo 
results from the exact values, as a multiple of the error in the mean. The second line shows 
the x 2 P er degree of freedom. Numbers in bold type indicate results which should occur 
with a statistical probability of less than 0.001. This table shows "bad" or "very bad" 
generators, grouped as to whether they failed the test at the level of 10 6 (very bad) or 10 7 
(bad) sweeps. 
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Table 2: As for Table 1, except here the number of sweeps is 5x10 for the Metropolis and 
Wolff algorithms, and 10 for Swendsen-Wang. This table shows "good" or "very good" 
generators, where the first (good) group of generators failed some tests at this level, while 
the second (very good) group passed all tests. 
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Table 3: Percentage deviation of the Wolff Monte Carlo results from the exact values for the 
energy and specific heat of the 2-d Ising model using different random number generators 
based on the lagged Fibonacci generator F(43,22,0). The binary operations tested were — , 
*, ©, and subtract-with-carry (SWC). A Weyl generator was also added to SWC (SWCW) 
and to F (43, 22,-) (Weyl). 
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Table 4: Percentage deviation of the Wolff Monte Carlo results from the exact values 
for the energy and specific heat of the 2-d Ising model using the standard 2-tap lagged 
Fibonacci generator F(55,24,0) and the 4-tap generator F(55, 24,16, 8,0). 
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Table 5: Deviation of the Wolff Monte Carlo results from the exact values, as a mul- 
tiple of the error in the mean, using the lagged Fibonacci generators F(250,103,©) and 
F(250,103,+). The 16 2 results are for 10 6 sweeps per run, and the 128 2 results are for 
8.5 x 10 4 sweeps per run. 
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Figure 1: Relative deviation in the Monte Carlo result for the specific heat of the 2-d Ising 
model, for the Wolff algorithm using a lagged Fibonacci generator. Each point denotes a 
different lag and a different binary operation for the random number generator. 
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